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^i. ' Abstract. We use a dynamical systems approach to study Bianchi type VIo 

^^ ' cosmological models containing two tilted 7-Iaw perfect fluids. The full state 

Z^ ^ space is 11-dimensional, but the existence of a monotonic function simplifies 

, the analysis considerably. We restrict attention to a particular, physically in- 

P3 ' teresting, invariant subspace and find all equilibrium points that are future 

^ ' stable in the full 11-dimensional state space; these are consequently local at- 

^~i tractors and serve as late-time asymptotes for an open set of tilted type VIq 

^■~^ , models containing two tilted fluids. We find that if one of the fluids has an 

f^ • equation of state parameter 7 < 6/5, the stiffest fluid will be dynamically 

' insignificant at late times. For the value 7 = 6/5 there is a 2-dimensional 

I ^ bifurcation set, and if both fluids are stiffer than 7 = 6/5 both fluids will 

►^ , have extreme tilt asymptotically. We investigate the case in which one fluid 

^"^ ' is extremely tilting in detail. We also consider the case with one stiff fluid 

*vj ' (7 = 2) close to the initial singularity, and find that the chaotic behaviour 

I ^ which occurs in general Bianchi models with 7 < 2 is suppressed. 

o: 
^, 

^ I 1. Introduction 

O . 
Q^. In [1] the asymptotic late-time behaviour of general tilted Bianchi type VIq 

s" I universes with a tilted 7-law perfect fluid was analysed. All of the future stable 

bi)' equilibrium points for various subclasses of tilted type VIq models, as well as for 

^ ' the general tilted type VIq models, were found. In particular, for the particular 

•jr^ . value of the equation of state parameter, 7 — 6/5, there exists a bifurcation line 

1^^ ' which signals a transition of stability between a non-tilted equilibrium point to an 

S . extremely tilted equilibrium point, confirming the observation in [2], which followed 

■ ■ ■ ' up on the earlier work on tilted fluids [3-9]. 

In this paper, we will generalize this work and use the dynamical systems ap- 
proach [10,11] to analyse the Bianchi type VIq model containing two (in general 
both) tilted (and separately conserved) 7-law perfect fluids [12]. We shall use the 
formalism employed in [1] , and pay particular attention to some of the cases which 
were not analyzed in detail there. In particular, we shall study the two physically 
important cases of when the second fluid is a comoving stiff fluid (i.e., F = 2 and 
F = in the notation below) , which represents a scalar field close to the singularity, 
and when the second fluid is of extreme tilt (i.e., V = 1). The study of extremely 
tilted fluids is partly motivated from the brane-world scenario where gravitational 
waves in the bulk can, through the interaction of the bulk Weyl tensor with the 
brane, give rise to extremely tilted fluids as seen from a brane point of view. 

The Bianchi type VIo model is not the most general model but is sufhciently 
general to account for many interesting phenomena. Bianchi type VIo models with 
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various sources of matter has been studied. The type VIo model with a comoving 
perfect fluid was considered in [13], with a magnetic and a 7-law perfect fluid 
in [14], and with a tilted 7-law perfect fluid in [1]. In all of these studies there 
has been evidence that the type VIq model is asymptotically self-similar at late 
times; the generic solution will approach an equilibrium point in the dynamical 
system. In the two tilted model the state space is 11-dimensional, and obtaining 
all of the equilibrium points is an extremely laborious task. However, guided by 
previous work we concentrate on the equilibrium points in a particularly interesting 
subspace; we will then show that these equilibrium points are stable in the full 11- 
dimensional state space. 

The paper is organised as follows: Next, in section 2, we write down the equa- 
tions of motion using the orthonormal frame formalism and discuss the various 
invariant subspaces. In sectional we discuss the (relevant) equilibrium points and 
their stability, and present some general results regarding the asymptotic behaviour 
of the tilted type VIq model. In section 0] we consider the case in which one fluid 
is extremely tilted. In sectional we consider the case of a stiff fluid. Close to the 
initial singularity a scalar field is essentially massless and can be modelled by a 
comoving perfect fluid with a stiff equation of state parameter 7 = 2. This will 
enable us to consider the initial singular regime in this case. Finally we summarise 
our results in sectional 

1.1. Preliminaries. In this paper we shall assume the existence of two (in general 
both tilting) perfect fluids. We shall denote the second fluid using tildes (e.g., f2, V), 
and use the notation F = 7, where by convention we assume that F > 7. We shall 
also assume that each fluid is separately conserved. 

The energy-momentum tensor for each tilted perfect fluid is 

(1-1) T^^ = {p + p)u^LU„+pg^„, 

where u^ = (cosh6',sinh6'c'') is the fluid velocity and 6 its rapidity. The spatial 
vector c° is chosen to be a unit vector in the tangent space of the surfaces of 
homogeneity; i.e. c'^Ca — 1. We will further assume that the each fluid obeys the 
barotropic equation of state, 

(1.2) p=(7-l)p. 

In terms of the unit normal vector u = Bq to the group orbits the energy-momentum 
tensor takes the imperfect fluid form, where p,p^qa,i^ab can be explicitly written 
down. For the Bianchi cosmologies we can always write the line-element in canonical 
form with respect to the structure constants C\^ of the simply transitive Bianchi 
group type under consideration, which depend only on time [15]. For the type 
VIq model, Ac = and Uab has two non-zero eigenvalues with opposite sign. This 
implies that we can choose a frame such that the structure constants can be written 
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ab = 0. 



Furthermore, the type VIq has n^ < r? . The equations of motion can now be 
written down in terms of expansion- normalised variables. 
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2. Equation of motion 
We define the expansion-normalised variables: 



(2.1) 
(2.2) 



-2i;+ \/3Si2 vSSis 

Sah = V32]i2 S+ + V32]_ V3S23 

\/3Si3 V0E23 S+ - voS- 

iV22 = iVag = V3iV, iV23 - N32 - \/3iV. 



We also introduce the three-velocity V by 



(2.3) 



sinh 9 = 



V 



VT^^v^ 



, o<y < 1. 



It is also convenient to introduce a parameter A defined by iV = AA'^, which makes 
it possible to solve for the function Ri (where Ra is the local angular velocity of 
a Fermi-propagated axis with respect to the triad Sa) and eliminate it completely 
from the equations of motion. For the Bianchi type VIq model this parameter will 
be bounded; i.e. A^ < 1. 

For the expansion-normalised variables the equations of motion are: 



(2.4) 
(2.5) 

(2.6) 

(2.7) 

(2.8) 

(2.9) 
(2.10) 



Si 



E'_ = 



y' — 



((Z-2)S++3(S22 + s23)-2iV2 

{q - 2)S_ + \/3(S?2 - S?3) - 2V3AS_E23 



S', 



13 

y' 

^23 

N' 
A' 



V^l^ I 2 2\ 

("2 - ^3) - 



2G+ ^ ^ ■'' 2G 
q - 2 - 3S+ - ^/3E_) S12 



^^^^{ii^~vl) 



^/3(AE_+S23)S 



13 



-V\V2 



V37 » 
q-2-3SH 



— ^j — W1W2 

G+ 
-V3E_')Ei3 



\/3(-AS_+S23)Si 



-W1W3 



-W1W3 



H — T^ — 1^ „ 

G+ G4 

(q - 2)^23 - 2\/3A7V2 + 2^/3AI]^ 



2^/3El2S 



12^13 



-W2U3 



W2W3 



V37SI 
' G+ '- ' G 

g + 2E+ + 2V3S23AI N 
2\/3E23(l - A') 
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The equations for the fluids are (similar equations for the second fluid): 

(2.11) n' = A|2g-(37-2) + [2g(7-l)-(2-7)-75]y2| 

v[ = (T + 2S+) vi - 2V3E13W3 - 2%/3Si2V2 

(2.12) -V3N{vl-vl) 

^2 - (T-E+- VSS-) V2-\/3(E23 + AE_)v3 

(2.13) +V3XNV1V3 + V3NviV2 

v'3 - (t-I]+ + V3I]_)v3-%/3(S23-AS_)w2 

(2.14) -V3XNviV2 - VSNviVs 

V(l~V^) 

(2.15) V = ,J(^„,);2 [(37-4)-5], 



where 



2 i + (7-i)y2 2 i + (r-i)y2 



"^2 ^^2 I y^2 I \^2 I v^2 I v^2 

„6 



S — E^ + El + ^2^2 + ^13 + ^23 



3' 



ta2 2 I 2 I 2 



G+ = l + (7-l)l^^ 
..... ^ (37-4)(l-y2) + (2„7)y2^ 
^'•^^) ^ - l-(7-l)V^2 ■ 

These variables are subject to the constraints 

(2.17) 1 = T.^ + N^ + n + n 



(2.18) = 2E_iV- 



(2.19) = -(I]i2 + Ei3A)iV + 

(2.20) = (Ei3 + E12A) iV + 



7rii;2 TD,V2 



G+ G4 
7rit;3 r51w3 



G+ G+ 

The parameter 7 will be assumed to be in the interval 76 (0,2] and we can without 
loss of generality assume that 7 < F. 

Eqs H2.18(l - H2.2()(l essentially deflne Vi, and so the governing equations are a 
twelve-dimensional set subject to a single constraint (eq. (|2.17|) '): hence the dy- 
namical system is effectively eleven-dimensional. 

2.1. The state space. The constraint H2.17|l implies that I]±, E12, E13, E23, N, 
il and (l are all bounded. Combining the equations for N and N we get the 
equations 

(2.21) {N±Ny = (q + 2J:+ ± 2V3^23) {N±N) . 
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Thus if the initial data have N"^ < N'^, then this will hold for all times^. Hence, 
for type VIq, N will be bounded as well. The invariant subspaces A^ ± iV = 
correspond to Bianchi type II universes. We also require < V, y < 1 for physical 
reasons. This implies the bounds 

I]2_ + S2_ + Ef^ + S?3 + S^g + iV2 < 1 

N'^ < TV^. 

We have, therefore, the 15 bounded variables I]±, S12, S13, S23J N, N, Vi, Vi, fi 
and Q. However, due to 4 constraints, these are not all independent. Hence, the 
state space can be considered a subspace of a compact region in M.^^. There are 
also some discrete symmetries which are intrinsic to the type VIq geometry [1]. 

2.2. Invariant subspaces. Some of the physically interesting invariant subspaces 
are as follows (we will also assume that the boundaries are included): 

(1) T{VIo): The fuU state space of tilted type VIq. 

(2) F{VIo): S- = 0, S12 = -Si3, vi = 0, z;2 = V3. 

(3) T2{VIq): a Bianchi tilted type VIq with two tilt-degrees of freedom; E12 = 
S23 = iV = (or S13 = Y.23^N = 0), V2 = 0. 

(4) Ti{VIo): Tilted Bianchi type VIq with one tilt degree of freedom; S12 = 
Si3 = 0, i;2 = i'3 = 0. 

(5) B{VIo): Non-tihed Bianchi type VIq; T.^ ^ T.^ ^ Y.13 ^ V ^ V = 0. 

(6) T^{II): The tilted type II boundary; TV = ±N . 

(7) B{II): Non-tilted type II; iV^ = TV^ and V = f = 0. 

(8) B{I): Bianchi type I universes; N ^ N ^V ^V = G. 

(9) dvB{I): The Bianchi type I vacuum boundary; N = N = V = n = V = 
fl = 0. 

2.3. A monotonia function. Fortunately, the above system of equations pos- 
sess a monotonic function with makes it possible to determine an important result 
regarding the asymptotic behaviour of the two-fluid model. We define the two 
functions: 

(2.22, ff="-^"''':-:'. ,9-"-^"'*'"" 



Then the function, 

(2.23) X 



1 + (7 -1)1^2' '- i + (r-i)F 



2 



I3n + (3n 

where x is bounded and satisfies — 1 < X < 1, satisfies the differential equation 

(2.24) ^'^3^r~7)(l-x'). 

Hence, assuming F > 7, then x' > and x is monotonically increasing. From this 
we can deduce that (for x^ 7^ 1) 

(1) lim^^_oo X = ^1; i-S-, at early times pfl -^ 0, 

(2) lim^^+00 X = +1; i-G-, at late times (3n -^ 0. 



^This can also be seen by using A and CQ. I2.1Ui : N'^ < N'^ -^ \^ < I. 
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At early times, we have that /3il -^ 0. li /3 ^ (i.e., not extreme tilt), it follows 
that rj — * and the first (less stiff) fluid is dynamically negligible at early times. 
Indeed, if neither fluids have extreme tilt asymptotically, then liniT-^-oo fi = 0, and 
liniT-^+oo f2 = 0. Hence, asymptotically the dynamical behaviour is that of a single 
fluid cosmology. It is therefore of interest to study the behaviour for when (at least) 
one of the fluids have extreme tilt and under what circumstances this can occur 
asymptotically. Note that in the case of a comoving stiff perfect fluid (/3 = 1 and 
r = 2), the future asymptotic behaviour is governed by the flrst fluid (and is given 
in [1]), and the behaviour close to the initial singularity need not be oscillatory. 

The behaviour of the monotonic function x combined with the earlier work [1] 
also gives us a hint of where to flnd the relevant equilibrium points. In [1] we 
noticed that all the locally future stable equilibrium points occurred in the invariant 
subspace F{VIo). The monotonic function x tells us that the late-time behaviour 
is almost entirely dictated by the less stiff fluid. Hence, one would expect that all 
the locally future stable equilibrium points for the two-fluid model also lie in this 
invariant subspace. Since the total state space is 11-dimensional, this observation 
simplifles our analysis considerably. In the following we shall therefore only consider 
equilibrium points in F{VIo), but the stability analysis will be performed in the 
full 11-dimensional state space. 

3. Late-time behaviour 

The system of equations possesses a wealth of equilibrium points. In Appendix 
IBl we have listed all of the physically interesting ones in the invariant subspace 
F{VIo). 

By inspection of the equilibrium points and their eigenvalues, we see that they 
all play a role in the late-time behaviour for various values of the equation of state 
parameters. Based on the analysis in the Appendix, the equilibrium points that are 
locally stable to the future are summarized in Tabled these represent the possible 
future asymptotic states of the models. It should be emphasised that some of these 
equilibrium points have zero eigenvalues so when we refer to stability we mean in 
the technical sense; e.g., there is 1-dimensional equilibrium set which is normally 
hyperbolic. 

There is another thing worth mentioning regarding the second fluid. If fJ — > , 
then the second fluid will decouple and become dynamically insigniflcant. For these 
cases the tilt velocity V of the equilibrium points indicates the tilt velocity close 
to the equilibrium points. Only away from (but arbitrary close to) f2 = is ^ a 
meaningful quantity. 

Note that for 7 = 6/5 there is a 2-dimensional bifurcation set and for 7 > 6/5 
both of the fluids are dynamically signiflcant at late times. If 7 < 6/5, on the 
other hand, the stiffest fluid becomes dynamically insigniflcant and the late time 
asymptotic behaviour of the model is entirely dominated by the less stiff fluid. 

4. An extreme tilted perfect fluid 
In this case we assume that the second fluid is of extreme tilt; i.e., 
(4.1) V^ = l. 

This case was not explicitly studied in [1]. 
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r 


Future Attractor 


< 7 < 2/3 


7 < r < 4/3 

r = 4/3 
4/3 <r 


It (I) 
Ie{I) 


2/3 < 7 < 6/5 


7<r<4/3-(37-2)/12 
r = 4/3-(37-2)/12 
4/3-(37-2)/12<r 


C{Vh) 

CEiVh) 


7 = 6/5 


7<r 


S{Vh) 


6/5 <7 


7<r 


S{VIo) 



Table 1 . The future stable equilibrium points for various equation 
of state parameters. 



By a careful analysis of the equilibrium points in the invariant set F{VIi^) we 
can prove the following: 

Theorem 4.1. When V'^ = 1, all of the equilibrium points of the resulting dynam- 
ical system, of physical interest (in the invariant set F{VIo)) necessarily have the 
property that (i) V'^ = (i.e., the first fluid is comoving), (ii) V'^ = 1 (i.e., the 
first fluid is also of extreme tilt), or (Hi) (wi,i'2,W3) = —uj{vi,V2,V3) (i.e., the two 
fluids are aligned). 

This makes the resulting analysis of the equilibrium points and their stability 
relatively straightforward (and in many cases the results can be inferred from the 
analysis of the single fluid case [1]). 

For the extremely tilted invariant subspace, V = I, the equation for fl simplifies 
to (with r < 2) 



(4.2) 



n' = {2q~2-s)n. 



The extremely tilted fluid loses its dependence on the parameter F, and dynamically 
it behaves similarly to an anisotropic radiation fluid with F = 4/3.^ 

Close to the initial singularity we typically have that g w 2 (separated by down- 
ward spikes in the chaotic case). Hence, since \S\ < 2E, we see that Cl is decaying 
towards zero as we approach the initial singularity. Thus for F < 2, the extremely 
tilted fluid will become negligible close to the initial singularity (in the sense that 

n~^o). 

At late times we have that /317 — > 0. But since (3 = 0, the future evolution is more 
complicated; we cannot conclude anything about the late-time behaviour from the 
monotonic function alone. However, some general results can be obtained from the 
existence of equilibrium points and their local stability and using Theorem 14. II 

From Table ^ we see that when the second fluid is extreme, then the future 
attractors are: Isil), for < 7 < 2/3 and 4/3 < F, CeIVIo) for 2/3 < 7 < 6/5 
and 4/3 - (37 - 2)/12 < F, and S(yio) for 6/5 < 7 and 7 < F. But both X£(/) 
and Ce{VIo) have il = (and fl ^ 0), so that for 7 < 6/5 the second (extreme) 
fluid is dynamically negligible at late times. Only when 7 > 6/5 can the second 



Compare ea. l4.2i with a non-tilted radiation fluid with anisotropic stresses for which the 



evolution equation takes the form Qr' 
normalised anisotropic stress tensor. 



{2<j - 2 - Sa5n"'')r2r where H"'' is the expansion- 
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fluid be non-negligible, and in this case the future attractor is £{VIo), where both 
fluids are of extreme tilt. 

In principle, we should study the special case in which both fluids are of ex- 
treme tilt separately. However, the main results follow directly from the arguments 
above. Both extreme tilting fluids will be negligible dynamically at early times for 
appropriate values of the parameters and the initial singularity will be oscillatory in 
general. And at late times the local future attractor will be £{VIo) for 6/5 < 7 < F. 

Thus we have the following: In the Bianchi type VIq model the extremely tilted 
fluid becom,es dynamically insignificant (Ti — > Q) in the future if the non-extrem,ely 
tilted fluid has 7 < |. Note that this includes, for example, the dust case (for which 
7=1). 

We have identified the extreme tilting fluid to be the second fluid. In the above 
we assumed that 7 < F. The case 7 > F is not obtained by simply interchanging 
the two fluids in the case described above (due to the asymmetry in the choice of 
the extreme fluid), and care should be taken in describing the asymptotic dynamics 
in this case. However, this can be done using Theorem 14 . 1 1 and the local stability 
of the equilibrium points (see below) . 

4.1. Null fluid; Relation to branes. In brane-world cosmology matter fields and 
gauge interactions are confined to a four-dimensional brane moving in a higher- 
dimensional "bulk" spacetime. The five-dimensional Weyl tensor in the generalized 
Randall-Sundrum-type models [16] is felt on the brane via its projection, f^i,. In 
general, in the 4-dimensional picture the conservation equations do not determine 
all of the independent components oiS^u on the brane. In these models the gravita- 
tional field can also propagate in the extra dimensions. For example, at sufficiently 
high energies particle interactions can produce 5D gravitons which are emitted into 
the bulk. In some applications these gravitational waves may be approximated as 
of type N [17]. 

Let us assume that the 5-dimensional bulk is algebraically special and of type 
N. Then there exists a null frame such that the 5-dimensional Weyl tensor is given 
by [18]: 

(4.3) Cafccd = 4CHii^{ar7lfc4m^j. 

Defining n'^ to be the normal vector on the brane, the non-local stresses on the 
brane can be written as 

(4.4) Sab = Cui, {laijn\n^) - m\{lcn% [4«n=) - m^,(4n^) 

where Eaiyf^ = S°'a — 0. Using the projection operator on the brane and defining 
r as the projection of the null vector i], onto the brane, we consider the physically 
interesting case in which rih = 0, so that 4 is null. The four-dimensional projected 
Weyl tensor on the brane can then be written 

(4.5) S^^ = -(^ e44. 

This is equivalent to the energy-momentum tensor of a null fluid, which is formally 
equivalent to the energy-momentum tensor of an extreme tilted perfect fluid. Using 
a covariant decomposition of £^,y, the non-local energy terms are given by: 

(4.6) U = e(4M'^)^ Q^ = e(4u")4> ^M^' = ^k>^^-)- 
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The equations on the brane now close and the dynamical behaviour can be analysed. 

We can investigate the effect this type N bulk may have on the cosmological 
evolution of the brane. Since dynamically F '^ 4/3 in the isotropic case, for 7 > 
r = 4/3 we need to reverse the arguments given in the previous subsection. First, 
we note that if 7 > 4/3, then the null fluid is not dynamically important at early 
times; that is, the effects of the projected Weyl tensor will not affect the dynamical 
behaviour close to the initial singularity'^. This supports the result that an isotropic 
singularity is a local stable past attractor, and hence constitute the most likely 
initial conditions in a classical brane- world [19]. Also, if 7 < 4/3, the extreme 
tilting fluid is dynamically negligible to the future. 

Therefore, the null brane fluid is not dynamically important asymptotically at 
early and late times for all values of 7 of physical import, supporting the phe- 
nomenological analysis presented in [17]. 

5. Stiff perfect fluid 

In this section we consider the case that the second fluid is a stiff comoving 
fluid. Close to the initial singularity a scalar fleld is essentially massless and can be 
modelled by a comoving perfect fluid with a stiff equation of state parameter F = 2. 
The stiff fluid was not explicitly studied in [1] . This will enable us to consider the 
initial singular regime as r — > — cx) in this case. Unlike for non-stiff fluids with 
F < 2, for which the initial singular regime possesses an oscillatory, and very likely 
a chaotic, behaviour [6, 10,20], in this case there exists a global past attractor [21]. 

In the Appendix we have proven 
Theorem lA. 21 For 7 < F = 2, and 11 > we have that 



(5.1) lim f = 0, lim (y.^ + 0)^1. 

The overbar denotes an appropriate mean value (see ea. l|A.3|l '). Hence, the stiff 
fluid will be asymptotically co-moving at early times. Note also that T,^ + il = 1 
implies fl = N = 0. This means that the second fluid 7 < F will be dynamically 
insigniflcant at early times. 

For the stiff fluid case, F = 2, there exist equilibrium points given by (usually 
referred to as Jacobs disc) : 



;2_ + sig + n = 1, 


A = 


1, 


S- = S12 = Si3 = 


= n = 


= 0. 



(5.2) V^V^N 

The equilibrium points can thus be considered as the solid disc S^ + E23 < 1 in 
the (S+, E23)-plane. The centre of the disc represents the FRW model with a stiff 
fluid and the unit circle correspond to vacuum Kasner solutions. 

By linearising the system of equations with respect to the Jacobs disc, we get 
the following eigenvalues from the shear, N, A and fi equations: 

Ai,2 = 0, A3 = -2V3i;23, ^4^5 = -3i;+ ± V3E23 

(5.3) Ae = 2(1 + S+ + V3S23), A7 = -4\/3S23, A8 = 3(2-7). 



■^In the brane scenario the Friedmann equation contains quadratic matter terms p^ (compared 
to the general-relativistic p dependence), so that close to the initial singularity we can effectively 
replace 7 — » 27. This result is consequently consistent with [19]. 



10 ALAN A. COLEY & SIGBJ0R,N HERVIK 




Figure 1. Jacobs Disc. Equilibrium points inside the triangular 
region are stable to the past. 

The remaining (independent) eigenvalues come from the velocities Vi which, accord- 
ing to Theorem IA.2I have Vi ^ 0. Hence, there exists a triangular region inside 
the Jacobs disc in which the equilibrium points are locally stable into the past 
(see Fig^). This indicates that solutions including a stiff fluid are not generically 
chaotic as we approach the initial singularity. From the existence of the monotonic 
function %, we know that the less stiff fluid has (3Q — > as we approach the initial 
singularity. From the ea. H5.1|l we see that fi > implies i7 ^ at early times; i.e. 
we expect that the fluid with 7 < 2 will decouple from the equations of motion and 
the initial singularity is almost entirely determined by the stiffest fluid. In this case 
the chaotic behaviour, conjectured to exist for F < 2, can be avoided. 

Therefore to conclude: // one of the fluids is stiff (T = 2), there exist local past 
attractors corresponding to non-tilted solutions with a stiff perfect fluid. Further- 
more, the less stiff fluid will be dynamically insignificant at early times. 

6. Discussion 

We have studied Bianchi type VIq universe models containing two tilted 7- 
law perfect fluids using dynamical systems techniques. The state space is 11- 
dimensional. We utilized the existence of a monotonic function, analysed the equi- 
librium points in a particular important invariant subspace, derived Theorem 4.1, 
and subsequently found equilibrium points that are future stable in the full 11- 
dimensional state space and are consequently local attractors; the solutions, which 
are summarized in Table 1, serve as late-time asymptotes for an open set of tilted 
type VIq models containing two tilted perfect fluids. 

We then studied the case in which one fluid is extremely tilting in more detail 
(this analysis, plus the analysis of the stiff fluid case, completed the analysis of [1]). 
We determined possible future behaviour and discussed when the extreme fluid is 
negligible asymptotically. 

We then applied these results to a problem in brane-world cosmology. The four- 
dimensional equations on the brane are modified by the effect of the projection of 
the five- dimensional Weyl tensor on the brane, £^1/, and when the bulk gravita- 
tional field can be approximated as type N gravitational waves close to the brane, 
this implies that the form of £^^ is that of a null fiuid, which is formally dynami- 
cally equivalent to a perfect fiuid with extreme tilt. Therefore, we can model the 
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dynamical effects of Ef^^u in brane cosmology by assuming that the second fluid is 
of extreme tilt. We found that the effects of E^tv are not dynamically important 
(at least in this class of models) in the asymptotic regime close to the singularity 
(as expected from the work of [17]) nor to the future at late times for physically 
important values for the fluid parameters. 

Finally, we studied the important case of a stiff fluid. We found that in this case 
there is a local past attractor (5.2), so that unlike previous analysis (for models with 
7 < 2), there are no chaotic oscillations as the initial singularity is approached. 
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Appendix A. Some simple proofs 
Lemma A.l. For < 7 < F < 2, we have 
(A.l) q < 2, and q = 2^Y? + n = l. 

Proof. We can use the constraint equation and write 

q^2-2N^ - f(^,V)n^ f{T,V)h 
where 
(K2^ 1 3(2-7) + (57-6)V-^ 

Note that for < 7 < 2, < F < 1, the function / obeys /(7, V) > 0, and equality 
holds if and only if 7 = 2, V = 0. From this we can see that ea. (jA.l|) follows 
automatically. D 

As the initial singularity is approached the solutions typically become oscillatory 
and chaotic. In order to control this oscillatory behaviour it is useful to introduce 
an appropriate mean value. A, defined by 

(A.3) A{t) = —^ I " A{t)dt, To > r. 

To - T Jt- 

This definition of a mean value is very useful since it basically measures the inte- 
grated effect of a certain quantity. The possible sharp spikes appearing in the cur- 
vature variables as the solution oscillates are therefore integrated out. As r — > —00, 
these spikes can therefore be ignored as an integrated effect. 

In the case of a stiff fluid, F = 2, the following theorem is useful: 

Theorem A. 2. For 7 < F = 2, and Vt > Q we have that 



lim y = 0, hm [Y? + n] =1 



Proof. For a stiff fluid the equations for V and fl are: 
V' = i2-S)V, 



(A.4) n' 



(q-2) + (2-S)- 



n. 
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As shown above, q is bounded hy q < 2. Therefore, there exists an e(T) > such 
that q = 2 — e. Integrating the latter equation in l)A.4|l (after dividing by tl and 
replacing (2 — S)V with V'), we obtain (r < tq) 



To 



To 



(A.5) InO =-2/ edr + ln(l + y-') 

Since ln(l + V"^) is bounded and e > 0, the right-hand side is bounded from above 
by In 2. The left-hand side, on the other hand, is bounded from below by lnf2(To) 
(which is negative). Hence, both sides must be bounded. More precisely, there 
must exist a. 6 > such that H. > 5 > for t < tq. Moreover, boundedness implies 
that J^° edr is bounded, and hence, 

1 r" 

e = / edr -^ for t -^ — cx), 



which leads to 9 —> 2 at early times. Using ea. (|A.ip we thus have linv^_oo ( S^ + fi 
1. 

Furthermore, we note that S obeys the bound \S\ < 2S, which for < 17 < 1 
implies that V is strictly monotonically increasing. Since V is bounded, either 
V ^ or S —^ 2. The latter case can only happen if S — > 1, which implies that 
n — > 0. However, this case is excluded since fl > S > ior t < t^. The theorem 
now follows. 

D 

Appendix B. Equilibrium points 

The system of equations possesses a wealth of equilibrium points. Here, we 
will describe some of the interesting ones which are important for the late-time 
behaviour. 

B.l. Non-tilted + Non-tilted. 

B.1.1. J{I), FRW:. T? = N = \ = Q = V = V = n = Q,n = l, q= i(37 - 2), 

0<7,r<2. 

Eigenvalues: 

Ai,2,3,4,5 = -^(2-7), A6,7-i(37-2), A8 = ~3(r-7), Ag.icii = (ST - 4). 

B.l. 2. C{VIo), Collins Vh: I]_ = Eia = S13 = E23 = A = 1/ = V^ = fi = 0, 

S+ = -i(37-2), N^ = -|(37-2)(2-7), f! = f (2-7), q = ^(37-2), 2/3 < 7 < 2, 

o<r<2. 

Eigenvalues: 
Ai,2 = -J(2-7)(l±v/^^^), A3.4 = -J(2-7)(l±y^^^), 
A5 = -^(2-7), A6,7 = ^(57-6), A8--3(r-7), 
Ag = (3r - 4) + i(37 - 2), Aio.ii = (3r - 4) - ^(37 - 2) 
B.2. Non-tilted -|- Intermediately tilted. 
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B.2.1. Jt{I), FRW with a tilted fluid: j:^=N = X = = V = 0, n = l, 
q = i(37-2). 0<7<2. 
For the second fluid: 



4 - V 

r = -, n = 0, ui = 0, V2 = U3 = -y=- 



Eigenvalues: 



Al,2,3,4,5 = "2(2-7), h,7^^{3l-2), As = (37 -4), A9,io,ii=0. 

B.2.2. Ct{VIq), Collins type VIq with a tilted fluid: E_ = S12 = S13 = S23 = A = 
F - 0, S+ = -i(37 - 2), N^ = ^(37 - 2)(2 - 7), ^i = |(2 - 7), g = i(37 - 2), 
2/3 < 7 < 2. 

For the second fluid: 

r = ^ - ^(37 - 2), ii - 0, 5i - 0, S2 = i^3 - ^- 



Eigenvalues: 



3G 

Ai,2, 3,4,5,6,7 as for C{VIo), As = — ;— (87 - 6), Ag = 0, 

4Gi 



(B.l) Vn = -^(37-2)(l±^l-8i|-^j 

B.3. Non-tilted + Extremely tilted. 

B.3.1. Ie{I), FRW with an extremely tilted fluid: j:^ = N = X^O^V = 0, 
f^ = l, g=i(37-2). 0<7<2. 
For the second fluid: 

< F < 2, Cl^O, wi = 0, W2 = W3 = ^. 

v2 

Eigenvalues: 

Al,2,3,4,5 = -^(2 - 7), A6,7 = ^(37 - 2), 

As -(37 -4), \, = -^i^-^, Aio,ii-0. 

The zero eigenvalues correspond to the fact that the equilibrium point is a member 
of a two-dimensional set of equilibrium points. 

B.3.2. CEiVIo), Collins type VIq with an extremely tilted fluid: S_ = S12 = S13 = 
S23 = A = V = 0, S+ = -i(37 - 2), N^ = ^(37 - 2)(2 - 7), « - 1(2 - 7), 
<7 = i(37-2), 2/3<7<2. 
For the second fluid: 

< F < 2, h^O, wi = 0, V2^V3^^- 

v2 



14 
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Eigenvalues: 



3(2 — r) 

Ai,2,3,4,5,6,7 as for C{VIo), As = — — (57-6), 



4r 



A9 = 



3r_4+-(37-2) 



3, , / /II7- 18 

' A.o,n = -2(37-2) (^li^^I-^ 



B.4. Intermediately tilted + Extremely tilted. 

B.4.1. S{VIo), a "f = 6/5 bifurcation: 

7=i 0<r<2, q=^, ^+ = -1' S_=S23 = 0, Ei3 = -Ei2 
5 5 5 

V 1 

Wl = Vl = 0, V2=V3 = -r=, < y < 1, V2 =V3 = -^, 



(B.2) 

n 

(B.3) 

where 

A 



[12-25(1 + A)7V2] (5 + y2) 
50(1 -y2) ^ 

1 



n 



Ei2 = — v/257V2(i_A)-6, Ar = 



1 + A 



257V2(-2v^2 ^ 3;^ ^_ 1) _ g(5^2 _^ ;l) 
25(1-1^2) ' 

V5±3(T/ + 1)(1-A)V^ 
5^/I 



A2 - 4A + 



- { [(A2 - 4A + 7)\^ + (1 + A)(A - 5)] ' - 18(1 - A)^} 

[(A2 - 202A + 121)1/ + A2 + 98A - ll] ^ - 432A(A + 5)(A + 1) 
^ ~ A2-202A+121 

3 I [(A2 - 2A - 39)F + A2 + 28A + 15] ^ - 36(3 + 8A + 23A2 + 2A3)| 

^ ~ A2 - 2A - 39 

(B.4) 

In addition we have to make sure that the energy densities are positive and that 
E12 is real. This puts bounds on N: 

12 



S?2>0 



N^ < 



N^ > 



N^ > 



25(1 + A)' 

6(5t/2 + 1) 
25(21/2 _^3A+1)' 
6 



(B.5) __ , .. -_^^^^_^^ 

The bifurcation S{VIq) is thus parameterized by A and V provided that the bounds 
(|R5|l hold. 

The region is bounded by the following lines: 

(1) fi = 0: Given by the line bifurcation C{VIq) for a single tilted fluid (see [1]). 

(2) n^O: 



, 12 



(3) V = 0: 



^2^9 + 42A 3A2±3(1-A)VA2 34A+1 ^ <,<,,_ ^^^ 
25(1 + A)(A2 + 14A+1) ' - - 
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(4) V=l: 

.12 ^6 

.^' '. ^ ^^■. 

Eigenvalues: Using a combination of numerical and analytical work the eigenvalues 

seem to be: 

A -n Rea \- ^ A - 6(5r^6) 

Al,2 — U, -n,e^A3^4,5^6J — ~7, A7 — — ^ , 

6 

-^8 + ^^9 = -^10 + All = —7, Re(A8,9,l0,ll) < 0. 

5 
B.5. Extremely tilted + Extremely tilted. 

B.5.1. £{VIq): Extremely tilted type VIq. V = V = 1, vi = vi = Q, V2 ^ v^^ ~ 
V2 = ""3 = 5\/2, S_ = E23 = ^ = 0, S+ = -f, S12 = -Si3 = ^, TV = 2^, 

rj + fi = ^,g = f. 

Eigenvalues: 



6(57 - 6) , _ 3 ^^ ^ . ^^ ^ _ 3 

5(2-7) 

h.7 = -l I 1 ± n/l4 + 5\/2 I , A8 = 0, A9 = -^^^^~^^ 



^— i^' .v,^-^(i..vi9),.„.-5(,,.^), 



[1 
[2: 
[3: 

[5: 

[6: 

[7: 
[»: 

[9 

[10: 

[11 
[12: 

[13; 

[14 

[is; 
[16; 

[17: 
[is; 

[19 

[20; 

[21 



5 V ^ 7 ' 5(2 -r) ' 

-^ 1 l±i\/ll + 4V2 
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